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Abstract: We use the variational method for the Sturm-Liouville eigenvalue problem 
to analytically calculate some properties of holographic superconductors with Gauss- 
Bonnet gravity in probe limit. By studying the holographic p-wave and s-wave su- 
perconductors in (3+l)-dimensional boundary field theories, it is found that near the 
critical temperature, the critical exponent of the condensation is 1/2 which is the uni- 
versal value in mean-field theory. We also find that when Gauss-Bonnet coefficients 
grow bigger the operators on the boundary field theory will be harder to condense. 
These are in good agreement with the numerical results. 
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1. Introduction 

The AdS/CFT correspondence 0, provides a theoretical method to understand 
strongly coupled field theories. Recently, it has been proposed that the AdS/CFT 
correspondence can also be used to describe superconductor phase transitions 0, ||. 
Since the high T c superconductors are shown to be in the strong coupling regime, 
one expects that the holographic method could give some insights into the pairing 
mechanism in the high T c superconductors. 

There have been lots of works studying various holographic superconductors, see 
review ^,0, in which some effects such as scalar field mass, external magnetic field, 
and back reaction etc., have been discussed. Among those works, some universality is 
discovered. For instance, the critical exponents for the condensation near the critical 
temperature is found to be 1/2 which is a universal value in mean-field theory ||. 
However, this value of 1/2 is mostly obtained by numerically solving the holographic 
systems and then fitting the data. The frequently asked question is how to calculate 
this 1/2 from analytical study. In ||, the authors used the variational method for the 
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Sturm-Liouville (S-L) eigenvalue problem 1 to analytically calculate some properties of 
the holographic superconductors in a (2+l)-dimensional boundary field theory. They 
managed to obtain the critical exponent 1/2 by perturbing the system near the critical 
temperature. Afterwards, the authors in |TU] extended this S-L method to study the 
holographic s, p, d-wave superconductors in 3-dimensional boundary field theory. They 
found that this 1/2 is a universal value for these models. 

In this paper, we extended the S-L method to study holographic p-wave and s- 
wave superconductors in a (4+l)-dimensional bulk with Gauss-Bonnet gravity. It is 
not trivial to extend this analytical studies to (4+l)-dimensional case, because the 
dimension of space-time will affect the boundary behavior of the fields which coupled 
to gravity. And moreover it would further have an impact on the conformal dimension of 
the dual operators. For instance, in (2+l)-dimensional boundary field theory the charge 
density p oc T 2 where T is the temperature; However, in (3+l)-dimensional boundary 
field theory the relation is p oc T 3 . This discrepancy would make the analysis a little 
different from the one in |§. We will see it explicitly in the following context. Besides, 
we also calculated the ratio between the critical temperature and the charge density 
for various Gauss-Bonnet couplings. The ratio was greatly consistent with the previous 



numerical results [11, |12|. We also found that the greater the Gauss-Bonnet coefficients 



were, the harder the condensations to form, which was also consistent with the previous 



conclusions (jT^, [14], |To| . Although the conformal dimensions of the dual operators on 
boundary theory are different from those in ||, there still exists a universal critical 
exponents 1/2 as in the mean- field theory. The ratios between the condensation values 
and the critical temperatures matched the numerical results up to the same order. 

The paper is organized as follows: In Sec.(fJ) we briefly review the basic facts of 
Ricci flat Gauss-Bonnet- AdS black holes. In Sec.(Q), we will analytically calculate 
the ratios between critical temperature and the charge density as well as the critical 
exponents of holographic p-wave superconductors. The same procedure will be followed 
in Sec.(||) by studying holographic s-wave superconductors. We will draw a conclusion 
in Sec.flSf). 



2. Ricci Flat Gauss-Bonnet- AdS Black Holes 

The action of Gauss-Bonnet gravity in a 5-dimensional space-time can be written as 
S = J d 5 x^(R+^ + ^(R 2 -AR^R^ + BT^R^j), {2.1) 



1 For short, we will call this method used by jj) to be "S-L method" in the following context. 
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where 12/L 2 is related to the cosmological constant and L is the radius of the AdS 
space-time. The quadratic curvature term in this action is the Gauss-Bonnet term 
with a the Gauss-Bonnet coefficient. The Ricci flat solution of this action is the 5- 



dimensional Gauss-Bonnet-AdS black hole [IS]. The metric is described by 



1 r 2 

ds 2 = -f(r)dt 2 + -rr^dr 2 + -r^idx 2 + dy 2 + dz 2 ) 
f(r) L 2 



(2.2) 



f(r) 



r 

2a 



1-1/1 



4a 

Z 2 



(1 



ML 2 



where M is the mass of the black hole. The horizon is located at r 
and the temperature of the black hole is 



T 



nL 2 ' 



(2.3) 



Here we should notice that in the asymptotic region with r — > oo, 

,2 



f(r) 



r 
2a 



1 



1 - 



4a 
I 2 



Therefore, one can define an effective radius L e q of the AdS spacetime by 

2a 



T 2 



la 
L 2 



(2.4) 



(2.5) 



We can see from this equation that one has to have a < L 2 /4 in order to have a 
well-defined vacuum for the gravity theory. The upper bound a = L 2 /4 is called 
Chern-Simons limit. In the AdS/CFT correspondence, this asymptotic AdS spacetime 
is dual to a conformal field theory living on the boundary r — > oo. The temperature 
of the black hole is just the one of the dual field theory. For simplicity of making 
integration and making comparison to the previous numerical results, we will assume 
< a < L 2 /5 in the following context. 



3. Holographic P-wave Superconductors 
3.1 Basic setup 

In order to study the holographic p-wave superconductors in the probe limit, we will 
introduce an SU(2) Yang-Mills action in to the bulk theory, the Lagrangian density is 

J? p = -- A F«F^\ (3.1) 
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where, F* u = d^A a v - d v A* + e ahc A\A c v is the Yang-Mills field strength, a,b,c= (1, 2, 3) 
are the indices of the generators of SU(2) algebra, /i, v = (t, r, x, y, z) are the labels of 
space-time with r being the radial coordinate of AdS. The A ^ are the components of 
the mixed-valued gauge fields A = A^^dx^, where r a are the SU(2) generators with 
commutation relation [r a , r b ] = e abc T c . e abc is the totally antisymmetric tensor with 
e 123 = +l. 

Following Refs. [0, [l8|, |20|, , we choose the ansatz of the gauge fields as 



A(r) = (f)(r)T 3 dt + i)(r)r l dx. (3.2) 

In this ansatz we regard the U(l) symmetry generated by r 3 as the U(l) subgroup of 
SU(2). We call this U(l) subgroup as U(l)3. The gauge boson with nonzero component 
ip{r) along x direction is charged under Af = <f)(r). According to AdS/CFT dictionary, 
<f)(r) is dual to the chemical potential in the boundary field theory while ip(r) is dual 
to the x component of some charged vector operator J. The condensation of ip(r) 
will spontaneously break the U(l)3 gauge symmetry and induce the phenomena of 
superconducting on the boundary field theory. 

The Yang-Mills equations with the above ansatz ( |3.2| ) are 

dl<j> + ld r <f> - ^fcf> = 0, 

(3.3) 

^ + ( l r + ^fW +^ = 0. 



In order to solve the equations ( |3.3|) , we need the boundary conditions for the fields 
ip(r) and <f)(r). At the black hole horizon, it is required that <p(r + ) = for the U(l) 
gauge field to have a finite norm, and ip{r + ) should be finite. Near the boundary of the 
bulk, we have 

(j>(r) ->• n- p/r 2 (3.4) 
i){r) V (0) +^ (2) /r 2 . (3.5) 

\i and p are dual to the chemical potential and charge density of the boundary CFT , 
i/j^ and tp^ are dual to the source and expectation value of the boundary operator 
respectively. We always set the source tp^ to zero, as we want to have a normalizable 
solution. 

3.2 Relations between the critical temperature and charge density 

In this subsection, we will adopt S-L method, which was pioneered by §|, to analytically 
calculate the ratios between the critical temperature T c and the charge density p of the 
holographic p-wave superconductors. 
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Define z — — , the EoMs (|3.3|) will become 



(3.6) 

r + (- + + = o. 

where a prime denotes the derivative with respect to z. In the following we will scale 
L — 1. When the temperature approaches the critical temperature i.e., T — > T c , the 
condensation approaches zero, viz. ip — > 0, from eq. ( |3.6| ) we can get 



- w 0. (3.7) 



With the boundary conditions ( |3.4|) , we can get that near the critical temperature, the 
electric field behaves like 

<p(z) « Ar+(1 - z 2 ). (3.8) 

where, A = p/r\.. Near the boundary, we will introduce a trial function F(z) into V^- 2 ) 
as in P 

V>l~o~^ (2) 4~(£>4*W (3-9) 

where, -F(O) = 1,F'(0) = 0. Therefore, using Q3.81 ) and (|3.9| ) the EoM of ip(z) reduces 
to 



p 3 + 4a(5^-3)-3 > /l + 4a( g *-l) 

(1 + 4a(^ 4 - 1) - y/l + Aa(z 4 - l))z 

16az 2 , v , 2 4a 2 (^ 2 -l) 2 

+ =F(z) + A . v '= F(z) = 0. 

l + 4a(z 4 -l)- v/l + 4a(z 4 -l) ( y/l + 4a(z* - 1) - l) 2 

(3.10) 

Multiplying both sides of the above equation with 

T{z) = -^(v/l + 4«(z 4 -l)-l), (3.11) 
2y/a 

We will convert the eq. (|3.10| ) to be 

d („, A 8^ 5 T-i/ \ 4a 2 (z 2 - l) 2 T(z) \ . 

— T(z)F (z) + ; v = F (2) + A 2 . v F{z) = 0. 

dz\ K J V V y/l + 4a(z*-l) V(v / l + 4«(2 4 -l)-l) 2 ; 



-0(2) P(«) 



(3.12) 
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From the Sturm-Liouville eigenvlaue problem [22], the minimum of eigenvalues of A 2 
can be obtained from the variation of the following functional 



\ 2 [F(z)} 



j l Q dz{T{z)F\zf + Q{z)F{zf) 



(3.13) 



ft dzP{z)F{zf 
The trial function F(z) can be assumed to be 

F(z) = l-az 2 , (3.14) 

which satisfies the boundary condition. So A 2 (z) can be explicitly written as 

s(a, a) 



X\z) 



t(a, a) 



(3.15) 



where, 



s(a, a) = 60(4a - 1) ((4a - l)a 2 + 4a) log 



t(a, a) 



(1 - 4a)a 



15 



(2a + v/a")' 

80v/a^l2a + a (20aa - 3a + 32^1 - 4aa - 48a - 8Vl - 4a + 8) J. 

. r . / -4a + \/l - Aa + 1 \ ,, 9 

a + l)a 128(a + l)log — a 3/2 

\ 2 — 8a / 

^1 — 4a)a 



+ (l6(3a + 4)a 2 + 8(3a + 2)a-a)log 

1 (2a + a/cFT 

-2v/a f (a (32a (Vl - 4a - 30 log 2 + 5) + 56^1 - 4a - 45) - 2y/l - 4a + 2) a 2 
+5a (4a (32^1 - 4a - 96 log 2 + ll) + 16^1 - 4a - 13) a 

+20a (a (20V1 - 4a - 48 log 2 + 6) + y/l-Aa - l) J (3.16) 

For different values of a, we list the minimum value of A 2 and the corresponding 
value of a in Table (|l|). Because T = r + /n and A = p/r^_, we can easily get that when 
T ~ T c , there is 



(3.17) 



where, 7 = l/(vrA^f n ). 

In Table (|ID, we also list the analytical values and numerical values [[12| of 7 for 
different a. The differences between the analytical and numerical values are within 
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Table 1: The minimum values of A 2 and the corresponding a for different a are listed in 
the first three columns. The analytical and numerical values of 7 and C are shown in the last 
four columns. The differences between the analytical and numerical values of 7 is about 1% 
for various a's. The order for Canalytical and Cnumerical is the same. 



a 


a 


\2 

min 


^analytical 


^numerical 


Canalytical 


Cnumerical 


0.0001 


0.685 


16.748 


0.199 


0.201 


326.858 


499.036 


0.1 


0.677 


20.604 


0.192 


0.194 


337.735 


509.796 


0.2 


0.653 


28.571 


0.182 


0.183 


349.614 


519.583 



1%. We can find that when a grows, the ratio 7 will decrease, which means the critical 
temperature will get lower compared to the charge density. This implies that the greater 
the Gauss-Bonnet coefficient is, the harder the condensation to form. This statement 



is consistent with the previous numerical results |TT|, [12], [13], [14 
3.3 Critical exponents and the condensation values 

From last subsection we know that when T — > T c , the condensation value of the dual 
operator (J*) is very small. Using formula ( |3.9| ) we can convert the eq.( |3.6|) to be 



Az 2 F 2 



(3.18) 



where, A = (J*) /?+• Because A is small, considering (|3.8|) we can expand (j)(z) in A 
as 



^~\(l-z 2 )+A X (z) 



(3.19) 



From the boundary condition we can get x(l) = 0, and for simplicity we set x'(l) = 
0. The eq.( [3.18|) now becomes 



„ X ' _ x z 2 F\l - z 2 
X — - A 7 

z f 



(3.20) 



When temperature is away but very close to the critical temperature, we expand 
(z) near z = 0, 



P ,2 



P 



P 



/i = -g"(l — Z ) = -^-(l — Z 



X(l - z 2 ) + A X (z) 



+ 



X(l-z 2 )+A{x(0) + x'(0)z + 



X"(0) , 



+ •••)• 

(3.21) 
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Comparing the coefficients of z 2 term in the above formula 2 we can get 



4=*-4*"(o). 

r%_ 2 



(3.22) 



From the eq. Q3.2U ) we know that 

-X'(z) Az 2 F 2 



x"(o) 



f 



■X'{z). 



2->0 



(3.23) 



2->0 



Because the second term in the r.h.s of the first equality approaches zero when z — > 0. 



In order to calculate the value of [^-M] 
both sides of eq. ( |3.20| ) we get 



fortunately we find that multiply 1/z to 



z->0 



1 zF 2 (l-z 2 ' 
z J 

Integrate both sides of the above equation and note that = 0, we obtain 3 

rX'(z). 



(3.24) 



x :\ zF 2 a-z 2 ) 

— A / az- 



Z->0 



/ 



(3.25) 



Therefore, from ( |3.22| ) one has 



P_ 

r 3 + " 2 
From the form of A and f(z) we can obtain 



, A, f 1 J zF 2 (\-z 2 ) 
A+ —A / cfe ^ 



(3.26) 



P _ 1 WIT 
Xri 2rl 



(3.27) 



where B = j dz(zF 2 {z){l - z 2 )/g{z)) and g{z) = f{z)/r 2 + . So we can deduce that 
where ( = \^6tt 3 / \/B. 



2 Here, the analysis is a little different from M in which they compared the coefficients of z term. 
The main difference comes from the near boundary behavior of 4>{z) which deeply originates from the 
dimension of space-time. 

3 Comparing the z term in ( 3.21 ), we know that x'(0) ^ when z 0. This is consistent with 
(p|). 
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The analytical values and numerical values [|T!| of £ for various a's can be found 
in Table. ([!]). They match up to the same order. Besides, in both analytical and 
numerical cases, the ratio ( increases when a grows bigger, which implies that the 
condensation values for the dual operators will also increase. From the AdS/CFT 
dictionary, this operator can be interpreted as the operator for the paring mechanism, 
the bigger vacuum expectation value of this operator will make the condensation harder. 
This statement is also in agreement with the previous arguments [O, 12, 13, 14 . . 



4. Holographic S-wave Superconductors 
4.1 Basic setup 

In order to study the holographic s-wave superconductors, we consider the Maxwell 
field couples to a complex scalar field in the bulk. The Lagrangian density is 0, 

J? s = -\f^ - |W - iA^f - m 2 |Vf (4.1) 

where A^ir) is the vector field and m 2 is the mass square of the scalar field. As usual, 
we can choose the gauge of vector field as = (<fi(r), 0, 0, 0) which is consistent with 
the ansatz of a real scalar field ip(r) = ip(r)*. In the following, we also scale the AdS 
radius to be L = 1. 

The EoMs for <p(r) and ip{r) are 

diet) + J d r <P - 2 -f<P 

Near the boundary, <p(r) and ip(r) behave like 

<Kr) (4-3) 

*tr) - ^ + ( 4 - 4 ) 

where, /i and p have the same interpretation as those in Sec.@. A± = 2± a/4 + m 2 L 2 s 
is the conformal dimension of the dual operator J in the boundary field theory. We 
will also set ip~ = and (J) = Therefore, we will unambiguously define A = A + 
from now on. On the horizon, 0(r) = because of the well definition of the norm. 



(4.2) 



2^ = 0. 
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4.2 Relations between the critical temperature and charge density 

Then, we can define z = r+/r, the EoMs Q4.2|) become 

,„ 1 ,/ 2^M. 



fz* 



0. 



(4.5) 



where a prime denotes the derivative with respect to z. 

Because when T — > T c , the condensation is very small, ip ~ 0. Therefore, the EoM 
of 4>(z) becomes 



0. 



(4.6) 



With the boundary condition ( f4.3|) , we can deduce that near the critical temperature 
<fi(z) ~ Ar + (1 — z 2 ) where A = p/r^_. 

As in the last section, we can introduce a trial function F(z) into the asymptotical 
behavior of il>(z), 



where F(0) = 1, F'(0) = 0. Now, the EoM for scalar field reduces to 
F"(z) +p(z)F'(z) + q(z)F(z) + X 2 w(z)F(z) = 0. 



where, 



p(z) 



q{z) 



(l - v/4 (z 4 - l) a + l) (2A - 3) + 4a ((2A - l)z 4 - 2A + 3) 



z(^4(z 4 -l)a-74(z 4 -l)a + l + l 

(^3 (z 4 -l)a + l - l) - 4 (z 4 - 3) aj A 

-4 (z 4 - 1) a + v/4 (z 4 - 1) a + 1 - 1 
2m 2 a 



(4.7) 



(4- 



(4.9) 



y/A (z 4 - 1) a + 1 - 1 
4(z 2 -l)V 



+ (A - 1)A 



10(2;) = — 

(V4(z 4 -l)a + l-l 

Multiply both sides of (fO|) with 



2 • 



(4.10) 
(4.11) 



T{z) 



z 



2A-3 



2v^ 



(v/l + 4a(z 4 -l)-l), 



(4.12) 
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we obtain the following equation 



where 



T(z)F'(z) - Q(z)F(z) + X 2 P(z)F(z) = 



Q(z) = -T{z)q{z) 
P(z) = T(z)w(z) 



(4.13) 



(4.14) 
(4.15) 



Following the procedure used in the above section, we know that the minimum value 
of A 2 can be obtained from variation of the following functional 



\ 2 [F{z)} 



j 1 Q dz{T{z)F\zf + Q{z)F{zf) 
Ji dzP{z)F{zf 



(4.16) 



It is easy to find that F(z) = 1—az 2 and in the calculation we will assume m 2 = — 3/L 2 ff 
in order to make comparison with the previous results [TT|. So \ 2 {z) can be explicitly 
written as 

s(a, a) 



\\z) 



t(a, a) 



(4.17) 



where, 
s(a, a) 



768a 2 



2v^^ (96V1 - 4aa + 404a - 75) a 2 
+48 (14^1 - 4aa - 24a - 5^1-40 + 5) a + 144 (2^1 - 4aa + a) ^ 
+3(4a - 1) ^25(4a - l)a 2 + 144a^ sinh" 1 ^2^-^- 



t(a, a) 



VI - 4a 

480a 3 / 2 



( (-144a 2 - 28a + l) a 2 - 40a(8a + l)a - 10a(20a + 



+-(a + l) 2 v^log (-4a + y/l-Aa + l) 
(a + 1) (16a(4a + 1) + a (48a 2 + 24a - 1)) log (y/1 - Aa^a) 



1 



128a 

960(a + l) 2 log(2-8a)a 2 



1920a 3 / 2 

-15(a + 1) (I6a(4a + 1) + a (48a 2 + 24a - l)) log (2a + 

+ (920a 2 + 90a - 4) a 2 + 10a(212a + 13)a + 960(a + l) 2 log(2)a 2 



2„ 2-r'/3\ \ nan/ „ i 1 \2_, „,2 



+40(24a + l)a - 960(a + l) 2 a 2 r'(-)/r(-) - 960(a + l) 2 7 Eu iera 



(4.19) 
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Table 2: For m 2 = -3/L 2 ff , the minimum values of A 2 and the corresponding a for various 
a's are listed in the first three columns. The analytical values from S-L approach are very 
similar to those obtained from matching methods |ll]] with differences smaller than 4%. The 
values of £'s are matching up to the same order. 



a 


a 


\2 

^min 


7s -L 


7matching 


Cs-L 


Cmatching 


0.0001 


0.722 


12.155 


0.210 


0.202 


7.705 


6.369 


0.1 


0.706 


14.290 


0.204 


0.197 


7.913 


6.884 


0.2 


0.662 


18.091 


0.196 


0.193 


8.042 


7.596 



in which, T(z) is the gamma function and 7Euier is the Euler's constant. Because 
T = r + /ir and A = p/r+, we can easily get that when T ~ T c , there is 

T c = 7 (4.20) 

where, 7 = 1/(7tAJ£). 

We list the analytical results in Table (|2[). We can find that the values of 7 ob- 
tained from the S-L approach are very close to those calculated in [ITJ by the matching 
methods 4 . The differences are smaller than 4% for various a's. We also find that 7 
decreases when a grows implying that the condensation is harder to form, which is 
consistent with the previous results Jll, 13 1. 

4.3 Critical exponents and the condensation values 

When T — > T c , from the equation of (f>(z), we get 

1 2A6z 2A ~ i F 2 
0" _ Ltf = ±^ 1_ (4.21) 

z 9 

where, A = (J) 2 / 'r 2A ,g{z) = /(j;)/r 2 . Because near the critical temperature T c , A is 
small. We can expand <fi(z) in A near z — 0, 

± = X(l-z 2 ) + A X (z) + --- , (4.22) 
r+ 

Comparing the coefficients in z 2 terms, we can get from the above formula that 

rl 2 



^ = A - ^ X "(0) (4.23) 



Substituting ( |4.22| ) into eq. ( |4.21| ), we can get EoM for x(z) as 



1 z 2A ~ 4 F 2 (l - z 2 ) 
X" ~ -X 1 = 2X- U j (4.24) 

z q 



'Consult Eq. (3.27) in Ref.Q- 
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Multiplying l/z to both sides and following the procedure in last section, we obtain 



X"(0) = 

Therefore, Q123| ) becomes 



2^0 



...3 



2A / dz 



1 ^2A-5 F 2 (1 _ ^ 



(4.25) 



A(l + ^B) 



(4.26) 



where B = 2 rf^ 2A - 5 F 2 (l - ,2 2 )/#. 

Finally, using the formula of ^4. we get that when T — > T c , the condensation of the 
operator J behaves like 



T 



(4.27) 



where, ( = y/Q/\/B. In Table(Q) we list the value of ( in S-L method and those 
in matching method [|11[] 5 . We find that these two kinds of values are matching up 
to the same order. The ratio ( increases when a grows bigger, which implies that the 
condensation will be harder to form. This is consistent with the previous results 11, 13|. 



5. Conclusions 

In this paper, we have used the S-L method to analytically calculate the holographic 
p-wave and s-wave superconductors with Gauss-Bonnet gravity in the probe limit. In 
particular, we calculate the ratios 7 between the critical temperature and the charge 
density on the boundary. The ratios are in great agreement with the previous numerical 
results or those obtained from the analytical matching methods. We find that the 
ratio will decrease when the Gauss-Bonnet couplings a increases, which implies that 
the bigger a is the harder the condensation to form. This statement is in agreement 
with the previous conclusions. Besides, we also calculate the ratio £ between the 
condensation value and the critical temperature. Although the values of ratio ( is 
not so well agreement with the numerical results or those from matching methods, 
they are consistent up to the same order. Moreover, we find that when a increases, 
the ratio ( will also increase. This again implies that the bigger a is the harder the 
condensation to form, which also agrees with the previous conclusions. The universal 
critical exponents 1/2 in mean-field theory was also analytically confirmed in both 
p-wave and s-wave case of Gauss-Bonnet gravity. This critical exponents will not be 

5 Refer to Eq. (3.28) in Ref.Q 



affected by the dimension of the space-time or the Gauss-Bonnet coefficients a. This 
value 1/2 is an exact property of the Lagrangian for the matter field of systems, and 
has nothing to do with the gravitational backgrounds . 
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